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Abstract 

In  Global  Navigation  Satellite  Systems  (GNSS),  the  on-board  time  has  to  be  modeled  and  predicted 
in  order  to  broadcast  the  time  parameters  to  final  users.  As  a  consequence,  the  time  prediction 
performance  of  the  on-board  clocks  has  to  be  characterized. 

In  order  to  estimate  the  time  uncertainty  of  the  on-board  oscillator,  a  linear  or  parabolic  fit 
is  performed  over  the  sequence  of  observed  time  difference  and  extrapolated  during  the  prediction 
period.  In  1998  the  Centre  National  d’Etudes  Spatiales  (CNES)  proposed  specifications  of  orbit 
determination  and  time  synchronization  for  GNSS-2.  The  needs  of  synchronization  were  specified 
as  the  maximum  error  of  the  time  difference  prediction  from  the  extrapolated  fit. 

Using  our  work  about  the  estimation  of  uncertainties  in  time  error  extrapolation,  we  have 
translated  these  time  domain  specifications  into  a  noise  level  limit  or  an  Allan  deviation  limit.  Of 
course,  these  limits  depend  on  the  main  type  of  noise  for  integration  time  of  about  1  day  and  on  the 
type  of  adjustment  which  is  performed  (linear  for  cesium  clocks  and  quadratic  for  other  oscillators). 

A  table  summarizing  these  limits  is  presented.  These  values  are  compared  to  experimental  results 
obtained  with  different  types  of  oscillators  (quartz,  rubidium,  and  cesium). 


INTRODUCTION 

In  Global  Navigation  Satellite  Systems,  the  on-board  time  has  to  be  modeled  and 
predicted  in  order  to  broadcast  the  time  parameters  to  final  users.  As  a  consequence, 
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the  time  prediction  performance  of  the  on-board  clocks  has  to  be  characterized. 

In  order  to  estimate  the  time  uncertainty  of  the  on-board  oscillator,  a  linear  or  parabolic 
fit  is  performed  over  the  sequence  of  observed  time  difference  and  extrapolated  during 
the  prediction  period,  i.  e.  when  the  satellite  is  out  of  visibility.  In  1998,  CNES  pro¬ 
posed  specifications  of  orbit  determination  and  time  synchronization  for  GNSS-2.  The 
needs  of  synchronization  were  specified  as  the  maximum  error  of  the  time  difference 
prediction  from  the  extrapolated  fit. 

The  question  is  then:  how  is  this  maximum  error  related  to  the  noise  levels  of  the 
clock?  Despite  several  papers  deal  with  this  question  [1,  2,  3,  4,  5],  a  new  approach 
was  chosen  here  because  we  are  not  only  interested  in  the  asymptotic  trend  of  this 
maximum  deviation,  but  also  in  its  evolution  close  to  the  interpolated  sequence. 

In  this  paper,  we  will  call  Time  Interval  Error  (TIE)  the  differences  between  the  ex¬ 
trapolated  parabola  and  the  real  time  deviation  x(t).  By  definition,  the  TIE  samples 
are  then  the  residuals  to  this  parabola.  However,  in  the  following,  we  will  limit  the 
use  of  the  word  “residuals”  to  the  differences  between  the  interpolated  parabola  and 
the  real  time  deviation  x(t). 

The  TIE  is  due  to  two  effects:  the  error  of  determination  of  the  parabolic  parameters 
and  the  error  due  the  noise  of  the  clock.  Obviously,  both  of  these  errors  may  be 
positive  or  negative,  and  the  ensemble  average  of  the  TIE  is  equal  to  zero.  Moreover, 
it  can  be  easily  shown  that  the  statistics  of  the  TIE  is  Gaussian.  Consequently,  we  only 
have  to  estimate  the  variance  of  the  TIE  in  order  to  completely  define  its  statistical 
characteristics. 

Moreover,  the  removal  of  a  quadratic  fit  from  the  time  deviation  sequence  cancels  out 
the  non-stationarity  problem  of  very  low  frequency  noises  [6,  7],  and  the  variance  of 
the  TIE  (i.e.  the  “true  variance”)  converges  for  all  types  of  noise  without  considering 
a  hypothetic  low  cut-off  frequency. 

In  order  to  determine  an  estimation  of  the  TIE,  we  will  already  redefine  the  interpo¬ 
lation  method.  Then,  we  will  compare  the  equations  giving  the  theoretical  estimates 
of  the  variance  of  the  TIE  to  simulations  and  to  real  data. 

ASSESSMENT  PRINCIPLE  OF  THE  INTERPOLATION  AND 
EXTRAPOLATION  ERRORS 

Interpolating  Functions  j 

Rather  than  performing  a  classical  linear  or  quadratic  fit,  we  used  the  Tchebytchev 
polynomials  [8,  7]: 

*(t)  =  Po$o(<)  +  H1$i(f)  +  T2$2(<)  +  e(t).  (1) 

where  e(t)  is  the  noise,  i.e.  the  purely  random  behavior  of  x(t),  {4>o(t),<I>i(<),<&2(t)}  are  the 
first  3  Tchebytchev  polynomials  and  the  parameters  {Pq,  P\,Pi}  have  the  same  dimension 
as  x(t),  i.e.  time. 
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Estimation  of  the  Residuals 

The  variance  of  the  residuals  a2  may  be  estimated  by  [8,  7]: 

ae  =  aX  ~  Jj  (aPo  +  VPi  +  °Pa)  ■  (2) 


Correlation  of  the  Samples 

Obviously,  the  long-term  behavior  of  the  TIE  depends  greatly  on  the  type  of  noise. 
The  autocorrelation  function  Rx(t)  of  the  x(t)  data  contains  the  information  about  this 
type  of  noise.  The  power  spectral  density  (PSD)  Sx(f)  is  the  Fourier  Transform  of  the 
autocorrelation  function  Rx(t). 

Taking  into  account  /;,  the  low  cut-off  frequency  and  //>,  the  high  cut-off  frequency, 
the  autocorrelation  function  may  be  written  [7]: 

Rx{tj  -  U)  =  f  Sx(f)  cos  -  ti)}  df.  (3) 

■hi 

We  will  use  the  power  law  model  of  Sx(f ): 

o 

SM)  =  E  *«/“•  (4) 

a~  —  4 


Calculation  Method  of  the  TIE 

By  hypothesis,  we  consider  that  the  TIE  is  the  difference  between  the  true  time  devi¬ 
ation  x(t)  at  time  t,  and  the  extrapolation  of  the  parabola  (previously  estimated  from 
to  to  tjv-i)  up  to  this  time  t  >  t^-i: 

TIE(t)  =x(t)-Po$o(t)-Pi$i(t)-P2$2(t)  and  t  >  tN_i.  (5) 

Thus,  the  quadratic  ensemble  average  of  TIE  may  be  estimated  by: 


<TIE2(t))  =  <*2W)  +  (A2)*o(0  +  (A2)®iW  +  (A2)®3(<) 

-2  ^ar(^)Po)  $o(t)  +  ^ar(J)A)  $i(t)  +  (*(f)A)  $2 (t) 

+2  [(AA)  +  (AA)  MhMh  +  (A  A)  ®i(*)M*) 


(6) 


Consequently,  for  each  type  of  noise,  we  have  to  know: 


•  (x2(t))  =  Rx(t ),  the  autocorrelation  function  of  x(t)  ; 

•  the  3  variances  ^A2)  =  aPi  > 

•  the  3  covariances  ^AA)  =  Cov(Pi,Pj):  actually,  only  (^PqP^  ^  0; 

•  the  3  covariances  ^ x(t)P, ^  —  Cov  (x(t),  Pj). 
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THEORETICAL  RESULTS 


Only  the  final  results  are  given  here,  but  the  calculation  details  are  available  upon, 
request  to  the  authors. 

In  the  following,  Tm  designates  the  measurement  time,  i.  e.  the  duration  of  the 
interpolated  sequence: 

Tm=Nr0.  (7) 

and  Tp  is  the  prediction  time  whose  origin  is  the  beginning  of  the  extrapolated  sequence 
(and  the  end  of  the  interpolated  sequence). 

Case  of  a  Quadratic  Interpolation 

Variance  of  the  residuals 

From  (2)  we  got  the  following  results  versus  Tm  =Nt0: 


White  FM  : 

2  /rp  \  37T 2k-2Tm 

<rATm)*  35  • 

(8) 

Flicker  FM: 

a2(T  )&  7r2^_3"i™ 

(9) 

Random  walk  FM: 

2/ti  \  ^k~^m 

<Xe(Tm)w  315  • 

(10) 

All  the  above  equations  were  obtained  under  the  assumption  N  »  1,  i.  e.  Tm  »  r0. 


Estimation  of  the  TIE  using  the  noise  levels 

The  theoretical  calculation  of  (6)  yields  the  following  variances  versus  the  measurement 
time  Tm  and  the  prediction  time  Tp: 


White  FM: 


(TIE2(Tm,Tp)) 


67r2fc_2Tm  (rX  _  T3 


35 


50^-  +  100-4-  +  6Mr  +  19^  +  1 

-t  rn  -i  m  -*  m  m 


(11) 


Flicker  FM: 


<TIE2(Tm,Tp)) 


7T2fe-3^ 

8 


rp&  rp5  rp 

192-^-  +  576-|-  +  692-2-  +  424-^-  +  136-f-  +  20-f  +  1 

fY\  I  rri.  2  m  I  m  -f  m  -t  771 


+96T3  111  VTm  +  Tp 


Lm 

/  rp4  rp3  rp2  rp  'S 

2^  +  7^  +  9^  +  5^-hl 

i  -Lm  -*■  m  1  w  , 


(12) 
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•  Random  walk  FM: 


<TIE2(Tm,Tp)) 


0_4t  rp3  (  y4  J-3  J-2  rp 

tK t+m4+mi+42i 


+  2 


(IS) 


All  the  above  equations  were  obtained  under  the  assumption  N  >  1.  Thus,  under  this 
assumption,  the  dependence  versus  N  cancels  out:  the  variance  of  the  residuals  only 
depends  on  the  length  of  the  interpolated  sequence  Tm,  and  the  variance  of  the  TIE 
only  depends  on  the  ratio  whatever  the  number  of  samples  N  is. 


Case  of  a  Linear  Interpolation 

Cesium  clocks  are  not  affected  by  a  quadratic  drift.  If  such  clocks  are  used,  we  may 
limit  the  fit  to  a  linear  interpolation.  In  this  case,  for  long  term,  the  main  contribution 
to  the  TIE  will  be  due  to  the  error  on  the  parameter  Pi  and  then  will  be  lower  than 
for  a  quadratic  interpolation.  It  is  then  strongly  recommended  to  extrapolate  the  time 
deviation  of  a  cesium  clock  from  a  linear  fit.  On  the  other  hand,  the  variance  of  the 
residuals  will  be  higher  because  a  quadratic  adjustment  remains  closer  to  the  time 
data  than  a  linear  one. 


Variance  of  the  residuals 


White  FM: 

_2 tm  \  ^ ,  ^2k-2Tm 
cre{Tm)  ~  j,- 

(14) 

Flicker  FM: 

a2(r  x  „  Tl 

VeVm)  «  g 

(15) 

Random  walk  FM1: 

<72(T  )  « 

aAm)  105 

(16) 

All  the  above  equations  were  obtained  under  the  assumption  N  »  1,  i.  e.  Tm  »  tq. 


Estimation  of  the  TIE  using  the  noise  levels 


•  White  FM: 


<TIE2(Tm,Tp))« 


47T  2k-2Tm 
15 


(17) 


•  Flicker  FM: 


<TIE2(Tm,Tp)) 


7T2fc_3r2 


rp4  rp3  rp2  rp 

12j£  +  24-4-  +  20^  +  8^  +  1 

-t  m  -I  m  m  TY), 


1  'Fhc  random  walk  FM  is  treated  here  for  homogeneity,  but  a  cesium  clock  is  never  affected  by  this  type  of  noise. 
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+2,n|1+f)  6i+6£+1 


rn3 

+2^  In 


t3  *“  Vrm  +  rp /  i 6t1  +  15rl  +  8 


(18) 


Random  walk  FM1: 


(TIE2(Tm,  Tp)>  *  ( 355.  +  39S.  +  11^  +  1 


105 


T3 

J-m 


T, £ 

m 


(19) 


All  the  above  equations  were  obtained  under  the  assumption  N  »  1,  i.e.  Tm  »  r0. 

Estimation  of  the  TIE  Using  the  Variance  of  the  Residuals 

The  relationships  (11)  to  (13)  and  (17)  to  (19)  needs  an  explicit  knowledge  of  the  noise 
levels  ka.  However,  for  very-long-term  interpolation  (several  days),  we  may  be  sui’e  of 
the  dominant  type  of  noise:  the  flicker  FM  for  a  cesium  clock  or  the  random  walk  FM 
for  a  quartz  oscillator. 

Thus,  the  variance  of  the  TIE  may  be  directly  estimated  from  the  variance  of  the 
residuals. 

As  an  example,  the  variance  of  the  TIE  for  a  random  walk  FM  and  a  quadratic  fit 
may  be  rewritten  from  (13)  and  (10): 

(TIE2(Tm,  Tp)>  m  2<r2  ( 450^.  +  690^- +  303^- +  42^  +  2  )  .  (20) 

\  1m  1m  Jm  1m  J 

However,  such  a  method  is  less  precise  than  the  use  of  a  correct  estimation  of  the  noise 
levels.  This  is  due  to  the  statistics  of  the  estimate  of  the  variance  of  the  residuals  (see 
section  “CONFIDENCE  INTERVALS”). 

Relationships  With  the  Time  Variance 

The  Time  Variance  (TVAR),  and  its  square  root,  the  Time  Deviation  (TDEV),  are 

commonly  used  for  time  analysis  [9].  Table  1  gives  the  relationships  between  the 

variance  of  the  residuals  and  TVAR(r)  for  an  integration  time  r  =  Tm  =  Nt0.  It  is 
interesting  to  notice  that,  for  a  given  type  of  noise,  the  ratio  cr^/TVAR(r)  is  constant. 


EXPERIMENTAL  VALIDATION 

Monte-Carlo  Simulations 

In  order  to  verify  the  equations  (11)  to  (13)  and  (17)  to  (19),  we  simulated  time 
deviation  sequences  of  different  types  of  noise  (white  FM,  flicker  FM,and  random  walk 
FM)  and  we  used  quadratic  and  linear  fits.  For  each  type  of  noise,  10,000  realizations 
were  calculated  with 

•  the  same  noise  level:  fe_2  =  1.4  •  10“4s,  &_3  =  3.3  •  10-8  or  k_4  =  5.0  •  10-12s-1  for  the 
quadratic  fit,  =  3.5  •  10_3s,  A:_3  =  4.8  •  10-7  or  fc_4  =  3.3  •  10_11s_1  for  the  linear  fit, 
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sAf ) 

k-2f~2 

k-3r3 

k-4f~4 

TVA  R  I'-r'l 

7 r&_2r 

7T2  [27  ln(3)  —  32  ln(2)]  ^_3r2 

ll7T4&_4T-3 

3 

6 

15 

Qua 

.dratic  interpolation 

Ji 

3ir2k-2T 

7r2fc_3T2 

7T4k-4T3 

35 

24 

315 

<?e 

9 

1 

1 

TVAR(r) 

35 

4  [27  ln(3)  —  32  ln(2)] 

231 

0.51 

0.18 

0.066 

TDEV(r) 

Li 

near  interpolation 

n-2 

27T2&_2r 

7T2k-3T2 

67T4&_4  7-3 

ae 

15 

9 

315 

2 

2 

2 

TVAR(r) 

5 

3  [27  ln(3)  —  32  ln(2)] 

77 

0"e 

0.63 

0.30 

0.16 

TDEV(r) 

Table  1:  Comparison  between  the  Time  Variance  TVAR(r)  and  the  variance  of  the  residuals  a2  for  different 
types  of  noise  and  for  quadratic  and  linear  interpolation. 


•  the  same  number  of  data:  65,536, 

•  the  same  number  of  data  taken  into  account  for  the  fit:  N  =  8640, 

•  the  same  time  of  estimation  of  the  TIE:  Tm+Tp  equals  to  {8640,  9900,  11350,  13000, 
14900,  17000,  19500,  22400,  25700,  29400,  33700,  38600,  44300,  50700,  58100,  65535}. 

The  noise  levels  were  chosen  for  getting  a  variance  of  the  residuals  equal  to  1. 

Figure  1  shows  the  curves  corresponding  to  the  square  root  of  equations  (11)  to  (13) 
and  (17)  to  (19)compared  to  the  standard  deviation  estimated  from  the  10,000  simu¬ 
lations.  The  simulations  exhibit  a  quite  good  agreement  with  the  theoretical  curves. 

The  asymptotic  behaviors  are  reached  at  the  end  of  the  log-log  graphs.  The  benefits 
of  the  linear  interpolation  are  obvious. 

Real  Clocks  (Quartz  USO  and  Atomic  Clocks) 

Figure  2  compares  the  long-term  behavior  of  a  real  ultra-stable  quartz  oscillator  (de¬ 
noted  “quartz  1”  in  table  2)  to  the  bounds  given  by  the  estimated  standard  deviation 
of  the  TIE  (the  square  root  of  Equations  (11)  to  (13)). 

The  data  from  the  oscillator  are  time  deviations  sampled  with  a  sampling  period  r0  =  10 
s,  obtained  at  CNES  Toulouse  (France)  with  their  own  reference  clocks  (Cs  HP  5071A 
option  001  and  H-maser  EFOS-16). 

Allan  variance  revealed  that  2  types  of  noise  must  be  taken  into  account:  white  FM 
(ho  =  10~22  s,  i.e.  2  =  2,5-  10-24  s)  and  random  walk  FM  (h- 2  =  1,5-  10-31  s-1,  i.e. 
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Clock 

ay(  24h) 

h—2 

(s-1) 

h- 1 

ho 

00 

Ce 

exp. 

(ns) 

Ce 

theo. 

(ns) 

&TIE 

exp. 

(ns) 

&TIE 
theo. 
(from  ha) 
(ns) 

cttie 
theo. 
(from  cre) 
(ns) 

Quartz  1 

1.8- 10"13 

0 

_ 2.2  ■  10-26 

7.5  •  10-23 

1.4 

1.4 

7.1 

6.2 

6.6 

Quartz  2 

2.9  -10"12 

1.4  10-29 

1.6  •  10-25 

0 

9.6 

9.2 

56.5 

52.0 

56.0 

Quartz  3 

3.0-  KT12 

1.4 -10-29 

6.4  •  lO-25 

0 

10.2 

11.0 

55.7 

59.0 

59.5 

Rb  1 

2.7-  1(T13 

1.2  - 10~31 

0 

5.3  - 10~22 

1.2 

1.3 

6.2 

5.7 

7.1 

Cs  1 

9.2  -  lO-14 

0 

0 

1.5  ■  10~21 

1.7 

1.6 

5.5 

5.5 

5.6 

Cs  2 

3.0  •  lO"14 

0 

2.1  •  10~28 

1.1  •  lO-22 

0.6 

0.5 

1.9 

1.6 

3.0 

Table  2:  Comparison  between  measured  and  estimated  values  of  cre  and  <ttie  for  several  clocks  with  a 
quadratic  interpolation. 


Clock 

°v(24h) 
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Table  3:  Comparison  between  measured  and  estimated  values  of  ae  and  cttie  for  two  cesium  clocks  with  a 
linear  interpolation. 


k- 4  =  4  •  10-33  s_1).  Thus,  the  bounds  of  Figures  2  were  obtained  by  using  the  square 
root  of  the  sum  of  (11)  and  (13).  This  fit  was  performed  over  a  24-hour  sequence,  but 
the  interpolated  sequence  was  shifted  along  the  90-hour  data  sequence.  At  each  step 
(to  =10s)  of  this  shift,  the  fit  was  extrapolated  over  3.5  hours.  The  TIE  measured  at 
this  instant  {Tp  =3.5  hours)  was  plotted  (solid  line)  and  the  interpolated  sequence  was 
shifted  again. 

The  TIE  bounds  of  the  left  figure  (dashed  lines)  were  estimated  from  the  noise  levels 
as  in  Figure  2.  In  the  right  figure,  the  TIE  bounds  were  estimated  from  the  variance  of 
the  residuals,  which  was  calculated  at  each  step  of  the  shift.  In  this  case,  we  assumed 
that  the  random  walk  FM  was  dominant  and  we  used  the  square  root  of  (20). 

The  experimental  TIE  curves  remain  generally  in  the  theoretical  bounds.  The  few 
moments  where  the  TIE  is  outside  the  bounds  is  fully  compatible  with  the  statistics 
of  TIE  (see  next  section). 

Table  2  shows  experimental  results  obtained  with  several  clocks  and  a  quadratic  in¬ 
terpolation.  The  noise  levels  were  estimated  from  Allan  variance  measurements.  This 
table  compares  “cre  exp.”  the  standard  deviation  of  the  residuals,  with  “<re  theo.J’ 
the  estimate  of  <re  obtained  from  the  square  root  of  (8),  (9),  (10)  and  the  noise  level 
coefficients. 

This  table  compares  also  “cttie  exp.”  the  measured  error  between  the  parabola  ex¬ 
trapolated  over  3.5  hours  and  the  real-time  deviation  of  the  clock  at  this  instant,  with 
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“<tt ie  theo.  (from  kaY  estimated  from  the  square  root  of  (11),  (12),  (13),  and  with 
“c ttie  theo.  (from  <re)”  estimated  from  the  square  root  of  the  expression  using  the 
variance  of  the  residuals  such  as  (20). 

Table  3  shows  the  same  type  of  results  but  limited  to  the  cesium  clocks  and  with  a 
linear  interpolation. 

Here  also,  the  agreement  between  measurements  and  estimates  is  quite  good. 

On  the  other  hand,  Table  2  and  3  confirm  that  the  standard  deviation  of  the  residuals 
is  higher  for  a  linear  interpolation  than  for  a  quadratic  one,  whereas  it  is  the  opposite 
for  the  standard  deviation  of  the  TIE. 

CONFIDENCE  INTERVALS  FOR  THE  INTERPOLATION  AND 
EXTRAPOLATION  ERRORS 

Figure  2  shows  that  the  real  TIE  may  be  outside  the  bounds  of  the  estimated  o-tib- 
Thus,  it  is  important  to  know  the  probability  for  the  real  TIE  to  be  inside  or  out¬ 
side  these  bounds.  Moreover,  it  may  be  useful  to  improve  the  TIE  bounds  by  using 
confidence  intervals.  This  may  be  achieved  by  studying  the  statistics  of  the  TIE. 

We  want  to  obtain  a  coefficient  cp  ensuring  the  following  confidence  interval: 

-cp  ■  (Ttie  <  TIE(Tm,Tp)  <  +<:p  •  <ttie  with  /?%  of  confidence .  (21) 

By  definition,  cp  follows  a  Student  distribution  with  v  degrees  of  freedom  [10].  Con¬ 
sequently,  for  building  the  confidence  interval  (21),  we  have  to  use  the  Student  coeffi¬ 
cients  for  cp.  These  coefficients  are  given  in  tables  [10]. 

However,  it  is  then  necessary  to  know  the  number  of  degrees  of  freedom  of  the  Stu¬ 
dent  distribution,  or,  and  this  is  equivalent,  the  number  of  degrees  of  freedom  of  the 
X2  distribution  followed  by  (TIE2(Tm,Tp)).  Obviously,  this  number  depends  on  how 
(TIE2(Tm,Tp))  is  estimated,  i.  e.  from  the  variance  of  the  residuals  or  from  the  noise 
levels  ka. 

Estimation  From  the  Variance  of  the  Residuals 

Equation  (20)  shows  that  (TIE2(Tm,Tp))  is  the  product  of  a  random  variable  aj  by  a 
constant  number.  Thus,  the  distribution  of  (TIE2(Tm,Tp))  is  the  same  as  the  distribution 
of  <Xg,  i.  e.  a  x2  distribution. 

From  Monte-Carlo  simulations,  we  observed  that  the  degrees  of  freedom  of  the  distri¬ 
bution  of  <t2  only  depend  on  the  type  of  noise  but,  curiously,  neither  on  the  number  of 
data  N,  nor  the  sampling  period  r0. 

We  measured  the  following  degrees  of  freedom  u  for  the  y2  distribution  of  <r 2  : 

•  White  FM:  v  »  8; 

•  Flicker  FM:  v  «  3; 

•  Random  walk  FM:  v  m  2. 
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1.1 

2.3 

k-3f~3 

3 

1.2 
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k- 4/"4 

2 

1.4 

4.3 

Table  4:  Degrees  of  freedom  of  and  (Ttie-,  and  confidence  coefficients  cp  for  70%  and  95%  versus  the  types 
of  noise  [10]. 


Therefore,  the  cp  Student  coefficient  must  be  chosen  with  these  degrees  of  freedom, 
according  to  the  type  of  noise.  Table  4  gives  these  coefficients  for  a  70%  confidence 
interval  (la-)  and  for  a  95%  confidence  interval  (2a-). 

Large  Estimation  R'om  the  Noise  Levels 

In  this  case,  the  degrees  of  freedom  of  (TIE2(Tm, Tp))  are  equal  to  the  ones  of  the  esti¬ 
mated  noise  level  %a  and  then  depend  on  the  accuracy  of  its  estimation.  For  instance, 
if  ka  was  estimated  by  using  the  Allan  variance,  its  degrees  of  freedom  depend  on 
the  length  of  the  sequence  and  on  the  number  of  sample  used  by  the  Allan  variance 
[11,  12].  If  this  sequence  is  long  enough,  the  degrees  of  freedom  may  be  much  greater 
than  the  values  obtained  from  the  variance  of  the  residuals.  The  degrees  of  freedom 
have  to  be  estimated  at  each  noise  level  measurement  [11,  12]. 

Therefore,  if  the  noise  levels  are  precisely  determined,  the  estimation  of  S-Tie  is  far 
better  by  using  this  method  than  from  the  variance  of  the  residuals. 


CONCLUSION 

The  first  application  of  this  work  may  be  the  selection  of  clocks  according  to  their  time 
stability  performances.  We  may  fix  a  limit  for  the  maximum  acceptable  ae  and  TIE 
for  a  given  interpolation  sequence  and  extrapolation  time.  Let  us  consider  that  for 
an  interpolated  period  Tm= 24h  and  for  an  extrapolated  time  Tp= 3.5h,  we  fix  :  <re  <2.1 
ns  and  TIE<5ns.  The  use  of  Equations  (8)  to  (13)  allows  us  to  translate  the  above 
specifications  into  specifications  on  the  noise  levels  ka  (for  instance,  for  a  random  walk 
FM,  these  specifications  become  &_4  <  3.7  •  10~33s_1).  Furthermore,  these  specifications 
may  be  translated  again  in  term  of  Allan  variance  over  1  day  (in  the  case  of  the 
random  walk  FM,  it  yields  <ry(24/i)  <  3  •  10-13).  Thus,  we  can  obtain  a  very  simple 
criterion  by  using  the  Allan  variance,  ensuring  that  the  specifications  for  cre  and  TIE 
will  be  respected. 

Besides  the  interest  of  this  method  for  navigation  satellite  systems,  it  may  be  used  for 
defining  a  new  method  for  very-long-term  stability  analysis. 

A  clock  may  be  continuously  measured  during  a  few  days  (e.g.  a  time  deviation 
measurement  with  a  sampling  period  of  1  minute  during  10  days).  From  these  data, 
the  noise  levels  of  this  clock  could  be  precisely  determined  and  a  quadratic  fit  could 
be  carried  out.  Thus,  if  the  clock  is  continuously  running  in  the  same  conditions,  it 
could  be  possible  to  extrapolate  the  difference  of  this  clock  with  the  parabolic  fit  after 
a  few  months  or  1  year. 


54 


This  analysis  could  be  helpful  to  low-accuracy  purposes  of  time  keeping,  for  instance 
for  industrialists  who  periodically  send  their  clock  to  an  accreditation  laboratory,  or 
for  applications  which  need  a  large  autonomy. 
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Figure  1:  Comparison  of  the  estimation  of  the  standard  deviation  of  the  TIE  calculated  from  the  equation 
(11),  (12),  (13),  (17),  (18),  (19)  (solid  lines)  and  estimated  over  10000  realizations  of  simulated  noise  for 
a  quadratic  interpolation  (o,  x,  and  □)  and  for  a  linear  interpolation  (*,  O  and  +).  In  order  to  use  the 
same  scale,  the  noise  levels  were  defined  in  such  a  way  that  the  variance  of  the  residuals  is  equal  to  one 
(&_ 2  =  1.4-  10_4s,  3  =  3.3  •  10~ 8  and  fc_4  =  5.0- 10~12s_1  for  the  quadratic  interpolation,  k- 2  =  3.5-10~3s, 

A;_3  =  4.8  ■  10-7  and  fc_4  =  3.3  •  10-11s-1  for  the  linear  interpolation).  The  error  bars  corresponding  to 
the  estimates  of  the  simulated  noises  are  too  small  to  be  plotted  on  this  graph.  The  right  plot  shows  the 
differences  between  the  theoretical  curves  and  the  simulation  points:  the  larger  error  is  equal  to  5%  but  most 
of  them  are  below  1%. 


Figure  2:  Estimation  of  the  TIE  for  an  ultra-stable  quartz  oscillator.  The  fit  was  performed  over  a  1-day 
sliding  window.  The  TIE  was  measured  3.5  hours  after  the  fitted  sequence  (solid  line).  The  TIE  bounds 
(dashed  lines)  were  estimated  from  the  noise  levels  (left)  or  from  the  variance  of  the  residuals  (right) . 
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